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In studying heavy-ion collisions at intermediate energies, it is important to know if ther- 
modynamical, i.e., thermal (momentum) and chemical (isospin), equilibrium can be reached 
in the heavy residues formed in these collisions. Although the early stage of heavy ion 
collisions must be described by microscopic models, a macroscopic treatment of the later 
stage in terms of temperature, volume and chemical potential becomes possible if thermo- 
dynamical equilibrium is established. Also, the interpretation of nuclear multifragmentation 
in heavy ion collisions either as a dynamical or as a statistical process depends on whether 
thermodynamical equilibrium is achieved in the collisions. Based on the assumption that 
thermodynamical equilibrium can be established in heavy-ion collisions, many statistical 
models have been developed, and they seem to be quite successful in describing the experi- 
mental data. A critical examination of thermodynamical equilibrium in heavy ion collisions 
is thus needed. A review on this subject can be found in Ref. |L|. 

In past few years, extensive theoretical efforts have been devoted to investigate the rate 
of thermalization of nuclear matter in intermediate energy heavy ion collisions. Using the 
Uehling-Uhlenbeck equation, momentum relaxation both in collisions of semi-infinite nuclear 
matter and in infinite nuclear matter have been studied • Collisions of finite nuclei have 
also been studied using the Boltzmann-Uehling-Uhlenbeck equation [^-0. On the other 
hand, only limited studies have been carried out for understanding the rate of chemical 
equilibrium in heavy ion collisions [||[10|- 111 both statistical and dynamical models it is 
usually assumed that chemical equilibrium is reached either instantaneously or as fast as 
thermal equilibrium. This is supported by deep inelastic heavy ion collisions at low energies, 
where the isospin degree of freedom has been found to first reach equilibrium [pTT| ^T3[| . Recent 



experimental studies |14}-|Ty| of the isotopic composition of intermediate mass fragments 
(IMF) and their angular distributions in heavy ion collisions have also shown that only at 
low energies (~ 30 MeV/nucleon) isospin equilibrium is reached before their emission. 

With recent advances in experiments using radioactive beams with large neutron or pro- 
ton excess, it has been possible to create a transient state of nuclear matter with appreciable 
isospin asymmetry. Knowledge on the rate of chemical and thermal equilibrium in such an 



asymmetric nuclear matter is thus needed [pl| . 



In this Brief Report, we shall report the results from a study of comparing the rate of 
chemical equilibrium to that of thermal equilibrium using an isospin-dependent Boltzmann- 
Uehling-Uhlenbeck (BUU) transport model. The model has been quite successful in pre- 
dicting and explaining a number of isospin-dependent phenomena in heavy-ion collisions at 
intermediate energies |) , 19[- 22 ] . The isospin dependence is introduced in the model through 



using different total and differential cross sections, Pauli blocking and mean fields for pro- 
tons and neutrons. A recent review of this model, especially its isospin-dependent aspects, 
is given in Ref . [[18] . 



In the model study of heavy ion collisions, a heavy residue can be identified as a collection 
of nucleons with densities higher than 1/10 of normal nuclear matter density. To characterize 
the degree of chemical equilibrium, we introduce the following quantity 

Mt) - (i) 

where (n/p) y>Q and (n/p) y<Q are, respectively, the neutron to proton ratio for positive and 
negative rapidity nucleons in the rest frame of the residue. If chemical equilibrium is estab- 
lished this quantity then has a value of one. We further define the isospin relaxation time 
Tj as the time when the quantity (A/(t) — l)/(A/(0) — 1) is 0.01, i.e., it is one percent from 
its equilibrium value. This time is an approximate measure of the rate at which the residue 
reaches chemical equilibrium. We note that this definition is somewhat different from that 
one would usually use, i.e., A(r/) = 1/e. 

For describing thermal equilibrium of the heavy residue, we use the quadrupole moment 
Q zz (t) in its rest frame, i.e., 

= J ^^Pl-pl~Pl)f(r,p,t), (2) 

where f(r,p, t) is the Wigner function from the BUU model calculations. Obviously Q zz = 
is a necessary, although not a sufficient, condition for thermal equilibrium. Similar to the 
definition of tj we define the momentum relaxation time t p as the time when Q zz (t)/Q zz (0) = 
1%. This quantity then measures the rate at which the residue reaches thermal equilibrium. 



Another important property of the heavy residue is possible existence of dynamical in- 
stability To study this phenomenon, we introduce the square of the adiabatic sound velocity 
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where (Ek) is the average kinetic energy per nucleon, a = -358.1 MeV, b = 304.8 MeV 
and a = 7/6 are the parameters corresponding to a soft nuclear equation of state. For 
v 2 < 0, a homogeneous nuclear matter is unstable against the growth of fluctuation, leading 
to dynamical instability or spinodal decomposition, 

We have studied the time dependence of (A/(i) — l)/(A/(0) — 1) and Q Z z{t)/Q Z z(0) com- 
pared to that of v 2 (t)/v 2 (0). They are shown in Fig. 1 for collisions of 40 Ca+ 124 Sn at an 
impact parameter of 1 fm and at beam energies of 25, 50, 150 and 300 MeV/nucleon. The 
system considered here has an initial A/(0) = 1.48. For collisions at a beam energy of 25 
MeV/nucleon the residue is found to be dynamically stable up to 300 fm/c during the col- 
lision. This time interval is long enough for both thermal and chemical equilibrium to be 
fully established as shown in the middle and lower windows of the first column. On the 
other hand, for collisions at beam energies above 50 MeV/nucleon a significant compression 
appears, and this is followed by expansion, leading into the adiabatic spinodal region after 
about 50 fm/c. We notice that at this time the heavy residue formed in the collision is still 
far from thermal and chemical equilibrium. Both the momentum and isospin asymmetries 
of the heavy residue are seen to change oscillatorily with time. The study on the nature of 
this oscillation in isospin asymmetry, in particular its possible connection to the giant dipole 
resonance in nuclei, is underway and will be reported elsewhere. 

We also notice from the middle window of the fourth column that nuclear transparency 
occurs at E/A=300 MeV. After spinodal decomposition the heavy residue quickly starts 
to break up into fragments and nucleons Although the isotopic contents of these 



fragments and nucleons depend on the emission angle, it is also strongly influenced by the 
neutron to proton ratio of the target and projectile in the entrance channel instead the 



average neutron to proton ratio of the combined system. This observation is consistent with 
recent experimental findings [P^I7|. 

Although chemical and thermal equilibrium are not completely established at beam en- 
ergies higher than the Fermi energy, it is still interesting to compare their relaxation times. 
This is shown in Fig. 2. In the left window, the comparison is made for 40 Ca+ 124 Sn colli- 
sions at beam energies from 25 to 300 MeV/nucleon and at an impact parameter of 1 fm. 
The momentum relaxation time is found to decrease with increasing beam energy. This is 
in qualitative agreement with that found in Refs. [§-|6|f|]. On the other hand, the isospin 
relaxation time decreases slowly with the beam energy. The shorter isospin relaxation time 
at incident energies below about 50 MeV/nucleon is in agreement with what was found in 
deep inelastic heavy ion collisions [|ITl -|T3|j. At higher incident energies the time for reaching 
momentum equilibrium is found shorter than that for isospin equilibrium. For example, a 
20 fm/c difference in the relaxation time is observed at E/A=300 MeV. In the right window, 
the isospin and momentum relaxation time as functions of the initial isospin asymmetry 
A/(0) = (^/p) P rojectiie/(^/p)target are compared for 40 Ca induced reactions on several isobaric 
targets of mass 124 at a beam energy of 300 MeV and at an impact parameter of 1 fm. It 
is seen that although the momentum relaxation is almost independent of the initial isospin 
asymmetry, the isospin relaxation time increases with the initial isospin asymmetry. 

The short relaxation time for isospin than momentum at low incident energies can be 
understood as follows. First, nucleon-nucleon collisions, which are responsible for momentum 
relaxation, are more likely to be suppressed due to Pauli blocking. Secondly, the repulsive 
symmetry potential for neutrons and the attractive symmetry potential for protons make 
preequilibrium emissions of neutrons more likely than protons in low energy collisions, as 
shown in Ref. P2 |, which thus enhances the isospin relaxation rate in the residue. On 



the other hand, in high energy collisions, Pauli blocking is less effective and the symmetry 
potential is also less important, leading thus to a shorter momentum relaxation time and a 
longer isospin relaxation time. 

We have also studied effects due to different forms of symmetry potential and the charge 



exchange reaction (jm — > np) on chemical and thermal equilibrium. They are found to have 
no discernible effects on both the momentum and isospin relaxation times. Only during the 
later stage of the collisions, they affect slightly the momentum and isospin distributions. 

In summary, using an isospin-dependent BUU model we have studied the isospin and 
momentum relaxation times in the heavy residues formed in heavy-ion collisions at inter- 
mediate energies. They are compared both with each other and also with the time for 
dynamical instability to appear in the heavy residue. We have found that chemical and 
thermal equilibrium can be completely established only at beam energies below the Fermi 
energy. At higher energies the dynamical instability sets in before either chemical or thermal 
equilibrium is achieved. Moreover, the isospin relaxation time is shorter (longer) than that 
for momentum at beam energies lower (higher) than the Fermi energy. 

This work was supported in part by the NSF Grant No. PHY-9509266 and the Robert 
A Welch Foundation under Grant A-1358. 
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FIGURE CAPTIONS 



Fig. 1 The development of dynamical instability (upper windows), thermal (middle 
windows) and chemical (lower windows) equilibrium in 40 Ca+ 124 Sn collisions at an 
impact parameter of 1 fm and beam energies of 25, 50, 150 and 300 MeV/nucleon, 
respectively. 

Fig. 2 Left window: isospin (open circles) and momentum (filled circles) relaxation 
times as functions of beam energy for the reactions shown in Fig. 1. Right window: 
relaxation times as functions of A/(0) of the reactions. 
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FIGURES 




FIG. 1. The development of dynamical instability (upper windows), thermal (middle windows) 
and chemical (lower windows) equilibrium in Ca+ 124 Sn collisions at an impact parameter of 1 
fm and beam energies of 25, 50, 150 and 300 MeV/nucleon, respectively. 
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FIG. 2. Left window: isospin (open circles) and momentum (filled circles) relaxation times 
functions of beam energy for the reactions shown in Fig. 1. Right window: relaxation times 
functions of A/(0) of the reactions. 
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